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Covering a graph or a lattice with non-overlapping dimers is a problem that has received consid¬ 
erable interest in areas such as discrete mathematics, statistical physics, chemistry and materials 
science. Yet, the problem of percolation on dimer-covered lattices has received little attention. In 
particular, percolation on lattices that are fully covered by non-overlapping dimers has not evidently 
been considered. Here, we propose a novel procedure for generating random dimer coverings of a 
given lattice. We then compute the bond percolation threshold on random and ordered coverings 
of the square and the triangular lattice, on the remaining bonds connecting the dimers. We obtain 
Pc = 0.367713(2) and Pc = 0.235340(1) for random coverings of the square and the triangular lat¬ 
tice, respectively. We observe that the percolation frustration induced as a result of dimer covering 
is larger in the low-coordination-number square lattice. There is also no relationship between the 
existence of long-range order in a covering of the square lattice, and its percolation threshold. In 
particular, an ordered covering of the square lattice, denoted by shifted covering in this work, has 
an unusually low percolation threshold, and is topologically identical to the triangular lattice. This 
is in contrast to the other ordered dimer coverings considered in this work, which have higher per¬ 
colation thresholds than the random covering. In the case of the triangular lattice, the percolation 
thresholds of the ordered and random coverings are very close, suggesting the lack of sensitivity of 
the percolation threshold to microscopic details of the covering in highly-coordinated networks. 


I. INTRODUCTION 

Percolation deals with the emergence of long-range 
connectivity in random systems, and is widely utilized 
in studying phenomena such as conductivity [1], mag¬ 
netism [2], porous fluid flow [3], permeation [4], jamming 
and glass transition [5] , gelation [6] and disease propaga¬ 
tion [7]. To that end, lattices and ordered networks are 
excellent model systems for studying percolation, with 
randomness introduced through using a known parame¬ 
terized probability distribution for occupying the bonds 
and/or sites of the corresponding network. The manifold 
in the parameter space of the probability distribution 
that corresponds to the emergence of long-range connec¬ 
tivity in the network is typically referred to as the perco¬ 
lation threshold, and has been computed analytically or 
numerically for a wide range of lattices and networks in 
different dimensions [8-33]. 

A particular class of lattices with unknown percolation 
characteristics are the dimer-covered lattices. In discrete 
mathematics, a dimer covering of a graph is a selection 
of edges of the graph so that each vertex is connected to 
exactly one such edge. Dimer coverings are of immense 
interest to fields such as discrete mathematics, statistical 
physics, chemistry and materials science. The problem 
of enumerating the number of distinct dimer coverings of 
a graph has received considerable interest in the litera¬ 
ture [34], especially in the asymptotic limit of infinitely 
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large lattices [35] . A conceptually related problem in sta¬ 
tistical mechanics is to determine the partition functions 
of on-lattice systems that can be mapped onto the dimer 
covering problem. For instance, there is a correspondence 
between the fully frustrated Ising model and the dimer 
covering of its dual lattice [36]. Dimer coverings of fi¬ 
nite honeycomb motifs are known as Kekule structures 
in chemistry, and are related to the number of ways that 
double bonds can be distributed on a polycyclic aromatic 
carbon backbone [37]. In materials science, dimer cover¬ 
ings arise in phases that assemble by dimers of a par¬ 
ticular building block. Such phases that are sometimes 
referred to as degenerate (guasi)crystals [38, 39] are ape¬ 
riodic and are the dimer coverings of the crystals [40] or 
the quasicrystals [41, 42] assembled by the corresponding 
monomers. These degenerate phases are stabilized by the 
pairing entropy, i.e., the entropy associated with distinct 
ways of covering a network with dimers. They are of par¬ 
ticular interest because of their superior mechanical prop¬ 
erties in comparison to their non-degenerate monomer- 
based counterparts [43, 44]. Degenerate phases of n-mers 
have also been observed in polymeric systems [45, 46], 
which can be conceptually considered as an extension of 
the dimer covering problem to situations where each ver¬ 
tex of the graph is part of a tree with n vertices. 

Despite the large body of work on enumerating the 
combinatorics of dimer-covered graphs and the partition 
functions of the associated model systems, very little is 
known about percolation on dimer-covered lattices. All 
existing studies have been carried out on lattices par¬ 
tially covered with dimers/Zc-mers by random sequential 
addition processes, and considering site percolation [47- 
55], or on systems such as interacting dimers [56], and no 
numerical estimate of the percolation threshold has been 
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obtained for fully dimer-covered lattices. Apart from the 
percolation threshold itself, a particular question of inter¬ 
est is the relationship between the percolation on a dimer- 
covered lattice, and the percolation on the associated un¬ 
covered lattice. As will be explained further in Section II, 
the constraint of dimers being non-overlapping will intro¬ 
duce some frustration in the percolation of the original 
lattice. In other words, the total number of bonds that 
need to be occupied for the lattice to percolate will in 
general be higher in the dimer-covered lattice. However, 
the extent of such frustration is not known a priori and 
can be different from lattice to lattice, and from covering 
to covering. 

Another point of interest is the inherent aperiodicity 
of random dimer coverings of a lattice. In other words, 
the topologically equivalent network of a random cover¬ 
ing (see Section II for details) is a random network itself. 
Therefore, the dimer covering process can be viewed as 
a novel way of generating random networks. Particu¬ 
larly, the construction of such networks is purely topo¬ 
logical and does not require the definition of a norm. 
Contrast this to other known examples of random net¬ 
works that need a notion of distance in their definition. 
Examples include the relative-neighborhood graph [30], 
Gabriel graph [19], Voronoi tessellation and Delaunay tri¬ 
angulation [17, 26]. Comparing the percolation charac¬ 
teristics of these covering-based networks with the norm- 
based random networks is also of interest. 

In this work, we study bond percolation on dimer cov¬ 
erings of the square and the triangular lattices using the 
Newman-Ziff algorithm [18]. We obtain random cover¬ 
ings of square and triangular lattice using a novel shuf¬ 
fling algorithm that is conceptually similar to an algo¬ 
rithms utilized for creating proton-disordered ice config¬ 
urations [57] . We assure the randomness of the generated 
lattices by computing appropriate correlation functions. 
We also compute the percolation threshold for several 
ordered dimer coverings of the square and the triangu¬ 
lar lattice for comparison. We refine our estimates from 
the Newman-Ziff algorithm using the epidemic growth 
method of Lorenz and Ziff [14]. We found thresholds to 
fairly high precision in order to distinguish the values for 
systems where Pc agreed to three or four digits. 

This paper is organized as follows. In Section II, the 
problem of bond percolation on dimer coverings of a lat¬ 
tice is outlined. The computational methods utilized for 
generating random coverings, and for measuring perco¬ 
lation thresholds are discussed in Section III. Section IV 
is reserved for results and discussions, and concluding 
remarks are presented in Section V. 


II. PROBLEM STATEMENT 

Studying percolation on dimer-covered lattices is a 
two-step process. Consider a lattice C that is comprised 
of Ng sites and bonds with Ng being an even inte¬ 
ger. The first step involves creating a dimer covering of 


£, denoted by Cd, in which (1/2)of those bonds 
are occupied permanently so that each site is connected 
to a single permanently occupied bond only. The second 
step involves defining the bond percolation problem on 
the arising network, La, by stochastically occupying the 
remaining Ni, — (1/2) Ag bonds with probability p, with 
Pc, the bond percolation threshold of Cd defined as [58]: 

Pc = sup{p : (/(p) = 0} (1) 

Here <p{p) is the probability of observing an infinite clus¬ 
ter in Cd- For a particular finite covering, the notion of 
an infinite cluster can be defined on a superlattice ob¬ 
tained from periodically replicating such a covering in all 
directions. This is the usual way of defining the perco¬ 
lation threshold as the occupation probability at which 
a percolating cluster emerges for the first time. Here, 
each cluster is comprised of sites that are connected to 
one another through a combination of permanently oc¬ 
cupied dimer bonds, and stochastically occupied regular 
bonds, and the definition of Pc is based on the bonds that 
are occupied stochastically. In other words, the percola¬ 
tion problem defined in (1) is for a new network in which 
each pair of sites connected through permanently occu¬ 
pied bonds is united into a single site. Throughout this 
work, we will refer to this new network as the topologi¬ 
cally equivalent network associated with the covering. In 
order to establish the relationship between percolation on 
the dimer-covered lattice, and on the original uncovered 
lattice, we define an equivalent percolation threshold as: 

Pc.eq = Pc (2) 

which is the percolation threshold of C conditioned that 
the permanently occupied dimer bonds in Ld are among 
the bonds that are chosen to be occupied in L. In 
other words, Pc,eq is the total number of bonds (perma¬ 
nently occupied dimer bonds, and stochastically occu¬ 
pied regular bonds) occupied when Ld percolates, with 
Ng/2Nb = 1/4 and 1/6 for the square and the triangu¬ 
lar lattice, respectively. Note that />con.d(p) < 'Pip) with 
4>con,di‘) and (/(•) the probability of observing an infi¬ 
nite cluster in L in the case of constrained and uncon¬ 
strained percolation, respectively. Therefore, Pc,eq > Pc 
with Pc the unconstrained bond percolation threshold of 
£. The ratio of k = Pc,eqlPc is a measure of percolation 
frustration due to the particular covering of C with non¬ 
overlapping dimers. This is a quantity of interest to us 
and will be computed and discussed in Section IV. 

It is necessary to emphasize that for any lattice £, a 
large number of distinct coverings are possible, and for 
any such covering, pc can be obtained from Eq. (1). We 
define the random covering percolation threshold as the 
ensemble average of Pc for all the coverings that can be 
obtained from £. 
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FIG. 1: Bond percolation in (a) parallel, (b) shifted and (c) staggered ordered coverings of the square lattice. The topologically 
equivalent lattices obtained by uniting the dimer sites are depicted in the bottom. Extrapolations to L —>■ oo have obtained 
using the linear fit of Pc vs. 


III. METHODS 


We compute Pc and Pc.eq for several coverings of the 
square and the triangular lattice, and we consider both 
ordered and random coverings. The particulars of the 
ordered coverings, as well as the procedure used for gen¬ 
erating random coverings are thoroughly discussed be¬ 
low. For three of the ordered coverings, Pc can be com¬ 
puted exactly through the triangle-triangle transforma¬ 
tion [21, 22], or through mapping onto the problem of 
inhomogeneous bond percolation on the square lattice. 
Even for these lattices, we perform the numerical esti¬ 
mation of Pc as a means of validating our computer pro¬ 
gram. In defining connectivity and covering, we apply 
periodic boundary conditions in both dimensions. In or¬ 
der to assess Pc in the limit of W —oo, four lattice 
sizes are considered: 128 x 128, 256 x 256, 512 x 512 and 
1024 X 1024. All random numbers are generated using 
the four-offset shift-register random number generator. 


i?(471,1586,6988,9689), described in Ref. [59]. 

A. Ordered Dimer-covered Lattices 

For the square lattice, three different ordered cover¬ 
ings are considered. In the parallel covering (Fig. la), 
individual dimer-covered rows are stacked on top of one 
another. The closely-related shifted covering (Fig. lb) 
is also comprised of stacks of dimer-covered rows. How¬ 
ever, the positions of dimers are shifted by one lattice 
site, every other row. The third covering is the staggered 
covering (Fig. Ic) and is comprised of horizontal and ver¬ 
tical stripes of the parallel covering. For the triangular 
lattice, the same types of coverings are considered. How¬ 
ever, here the parallel and shifted coverings correspond 
to the same lattice (Fig. 2). We will therefore refer to 
these two equivalent coverings as the parallel covering 
from this point onward. 

Among these five ordered coverings, the percolation 
















































































































































































4 



FIG. 2: Bond percolation in (a) parallel and (b) staggered ordered coverings of the triangular lattice. The topologically 
equivalent lattices obtained by uniting the dimer sites are depicted in the bottom. Extrapolations to L —>■ oo have obtained 
using the linear fit of Pc vs. 


A 



FIG. 3: The unit super-triangle of the parallel triangular cov¬ 
ering of Fig. 2a. 

thresholds of three can be computed exactly. As can 
be noted in the bottom of Fig. lb, the shifted covering 
of the square lattice is topologically identical to the tri¬ 
angular lattice, with its percolation threshold of Pc = 
2 sin ^ « 0.347296 • • • known analytically [8]. The par¬ 
allel covering of the triangular lattice is also topologically 
identical to a lattice that can be partitioned into super¬ 
triangles. Fig. 3 depicts one such super-triangle. There¬ 
fore, its percolation threshold can be computed from the 
triangle-triangle transformation [21, 22] i.e., by solving 


the equation P{ABC) = P{ABC). Here, P{ABC) is 
the probability that A, B and C in Fig. 3 are all con¬ 
nected, while P{ABC) is the probability that none of 
A, B and C are connected to one another. For the super¬ 
triangle depicted in Fig. 3, P{ABC) = P{ABC) yields 
the following algebraic equation: 

Ip^q^ + 9p^q^ + 5/(7 + p^ - q^ = 0 (3) 

with q = 1 — p. By solving (3) numerically, we obtain 
a percolation threshold of Pc ~ 0.237418 • • •. This also 
follows from the anisotropic triangular lattice threshold 
splitting one bond into three [8]. 

The third covering with a Pc known exactly is the 
parallel covering of the square lattice, which is topolog¬ 
ically identical to a square lattice with double bonds in 
one direction (Fig. la). Bond percolation on such a lat¬ 
tice can be obtained by mapping it onto the problem of 
inhomogeneous percolation on the square lattice. Sup¬ 
pose a square lattice in which the horizontal and vertical 
bonds have different occupation probabilities pi and p 2 - 
The percolation manifold of such a lattice is given by 
Pi + P 2 = 1- For the lattice of Fig. la, let pi = p and 
P 2 = 1 — (1 — = 2p — p^. The percolation thresh¬ 

old is therefore the root of — 3p -|- 1 = 0 which is 
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FIG. 4: A representative shuffling move of the algorithm used for generating random coverings of a lattice. Starting with an 
arbitrary initial covering, such as the parallel covering of (a), a site, such as point 1 in (b), is chosen randomly, and a closed 
path originating from 1 is constructed by carrying out a random walk, (c) The new covering is obtained by flipping the identity 
of permanently occupied bonds along the enclosed portion of the path (i.e., the portion starting and ending at 3). 


p= (3-v^)/2« 0.381966- 


B. Generation of Random Dimer-covered Lattices 

The algorithm that is used for generating random cov¬ 
erings of a lattice is based in creating closed paths of 
dimers, and shifting the permanently occupied bonds 
along that path. Starting with a given initial dimer cov¬ 
ering, such as the parallel covering depicted in Fig. 4a, 
a sequence of shuffling moves are performed. In each 
move, a closed path P is constructed as follows. First 
a site is randomly selected on the lattice, such as point 
1 in Fig. 4b. Both the selected site and its dimer pair 
(i.e., point 2 in Fig. 4b) are added to P. The next step 
is to choose a second dimer by randomly selecting one 
of the z — 1 unoccupied bonds ending in 2, and adding 
it to P. Here z is the coordination number of the orig¬ 
inal lattice. In Fig. 4b, for instance, 3 and 4 are added 
to P. This procedure is continued until P crosses itself 
at some point, e.g. by returning to 3 as in Fig. 4b. If 
the loop portion of P contains an even number of sites 
(which is alway the case for the square lattice but not the 
triangular lattice), the move is completed by flipping the 
bond identities along the closed portion of P, otherwise 
the move is rejected. For instance, the move depicted 
in Fig. 4b gives rise to the covering depicted in Fig. 4c. 
Upon performing a sufflciently large number of moves, 
the covering becomes fully randomized, by losing corre¬ 
lation with the initial covering, and eventually transforms 
to a covering such as the ones depicted in Fig. 5. The 
random coverings used in pc calculations are gathered ev¬ 
ery 1,000 steps, with each step consisting of Ng moves. 
For every lattice size/type, we gather a minimum of 100 
independent random dimer coverings. 

In order for our pc estimates to be reliable, it is nec¬ 
essary to establish the following two criteria: (a) In each 
random covering, no long-range spatial correlation should 


exist between the orientations of dimers, (b) The suc¬ 
cessive coverings generated by the algorithm are to be 
uncorrelated. In order to check for the first criterion, we 
compute a spatial orientational correlation function as 
follows. First, we assign to every site i an orientation Oi. 
In the square lattice, the sites participating in horizontal 
and vertical dimers are assigned the values of Oi = -bl 
and Oi = —1, respectively. In the triangular lattice, val¬ 
ues of -bl, —1 and 0 are assigned to horizontal, leftward 
and rightward dimers, respectively. Note that in a fully 
random covering, (oi) = 0 since it must be equally likely 
for a site to belong to dimers of different orientations. 
The spatial orientational correlation function is then de¬ 
fined as: 

9o{x,y) = {OijOi+x,j+y)i,j (4) 

with i,j,x and y the two-dimensional (integer) coordi¬ 
nates of the site on the lattice. In a fully-random cover¬ 
ing of a lattice, go{x,y) must decay exponentially along 
all arbitrary directions in the (x, y) plane. Fig. 6 depicts 
go{x,y) for a representative random covering of a 128- 
by-128 square (Fig. 6c) and triangular (Fig. 6c) lattice. 
Note the rapid decay of go{x, y) in all directions. go{x, y) 
is virtually zero beyond the third nearest neighbor shell. 
In contrast, the ordered lattices depicted in Figs. 1 and 2 
tend to have non-decaying spatial correlations in all di¬ 
rections (Fig. 6a-b). This shows that the algorithm uti¬ 
lized here creates coverings that are random, and that 
lack no spatial correlations in the orientations of dimers. 

To test (b), we compute the site orientational autocor¬ 
relation function given by: 

fs{t) = (Oj(T)Oi(T -b t))j,r (5) 

Here t corresponds to the number of steps, with each 
step comprising of Ng trial walks. We compute fg{t) for 
coverings of a 256 x 256 square and triangular lattice gen¬ 
erated by the above mentioned algorithm. The decay in 
fg(t) is very rapid and fg{t) becomes virtually zero for 
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FIG. 5: Schematic representation of random dimer coverings of (a) a square, and (b) a triangular lattice. The finite-size scaling 
of the computed percolation thresholds are given in the bottom. 


t > 2. This is much smaller than the 1,000-step sepa¬ 
ration used for gathering random coverings. Therefore, 
the distinct random coverings used in this work are truly 
independent and lack no correlation. 


C. Percolation Thresholds 

All percolation thresholds are computed using the 
Newman-Ziff algorithm [18]. In each trial, a sequence of 
randomly selected unoccupied bonds are occupied one at 
a time. Upon the addition of each bond, the largest clus¬ 
ter of connected sites is updated. The procedure is con¬ 
tinued until the largest cluster wraps around the lattice 
in either of the x oi y directions. The percolation thresh¬ 
old is then estimated by Pc = {Nm) /A^aii with being 
the number of bonds that have been occupied through¬ 
out the trial, and A^aii = Ni, — {l/2)Ns the total number 
of (initially) unoccupied bonds. The averaging is taken 
over different trials. For every lattice, whether ordered 
or randomly covered, we perform 10® trials, which gives 
a precision of three significant digits for p^- The random 
covering percolation thresholds are obtained through an 
additional level of averaging over a minimum of 100 inde¬ 
pendent random coverings. This increases the precision 
in the computed significant digits. 

The procedure outlined above gives Pc on a finite lat¬ 
tice. However, percolation transition is only precisely 
defined in an infinite system. There is a large body 
of work in the literature on how different estimates of 
Pc on finite lattices scale with the characteristic length 
scale of the lattice. A detailed overview can be found 
in Ref. [60]. It follows from general scaling arguments 


that \pc{L) — Pc,ao\ ~ with v = 4/3 for percola¬ 

tion in two dimensions. However, faster convergence can 
become possible depending on the system shape, bound¬ 
ary conditions, the lattice, and the numerical method 
that is used for estimating pc- In this work, we use both 

\Pc(L) -Pcooj ~ and \pc{L) -Pc,oo| ~ 

for hnite-size scaling. The latter scaling applies for a sys¬ 
tem with special symmetry, such as a square or rhombic 
shape, while the former applies for a general system. 

1. High-precision pc’s from Epidemic Growth 

To further increase the accuracy of the estimated per¬ 
colation thresholds, we carry out simulations using the 
epidemic growth method of Lorenz and Ziff [14] . In this 
method, a site v is randomly chosen on the lattice, and 
a cluster is grown from v via a Leath-type algorithm in 
which bonds are added at its periphery. Here, each bond 
is occupied with probability p and the process continues 
until there are no more growing bonds at the periphery, 
or when 2^"^ = 131,072 sites are already wetted by the 
growing cluster. Cluster size distribution data are ac¬ 
cumulated in logarithmic bins (2", 2"+^ — 1), which are 
then used for estimating P>s, the probability that the 
grown cluster is larger than s = 2". We perform this on 
sufficiently large lattice (8,192 x 8,192) so that the grow¬ 
ing clusters never wrap around the periodic boundary. 
Therefore, our estimates of P>s are not affected by the 
finite size of the system. According to the scaling theory, 
P>a scales as s^~'^ where r = 187/91 in two dimensions, 
and the finite-size scaling for small clusters is of the form 

s^-‘^P>s ~ A -h Bs-^ (6) 
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FIG. 6: Spatial correlation on: (a) parallel and (b) staggered, (c) random square and (d) random triangular dimer-covered 
lattices. 


with f2 = 72/91 [61, 62]. This relationship holds for 
sufficiently large When p is away from p^, there 

will be large deviations from the above behavior for large 
s. Fig. 7 depicts s'^~^P>s vs. for bond percolation 
on the staggered triangular lattice for values of p close 
to where the large-s behavior is linear. From this data 
we estimate Pc = 0.237497(2). This is consistent with the 
estimate of Pc = 0.237(1) obtained from the Newman-Ziff 
algorithm. 

For refining the error bars of the Pc’s estimated for ran¬ 
dom dimer coverings, we combine the above algorithm 
with partial rearrangements of dimers using a shuffling 
algorithm similar to the one explained in Section IIIB. 
For the square lattice, we continue the walks until the 
origin is revisited for the first time, and flip every bond 
along the path. For the triangular lattice, we flip bonds 
on a single intercepting loop, as long as the loop contains 
an even number of steps. Clusters are grown from ran¬ 
dom sites on the lattice. By having a sufficiently large 
lattice, we are able to grow an ensemble of uncorrelated 
clusters in different parts of the lattice without the need 
for completely rearranging the covering. We also use a 
technique of labeling occupied sites by the index of the 
number of runs, with the sites with indices less that that 
considered vacant. This way, we can avoid clearing the 
lattice between successive clustering events. 

For each covering, we start with a trial range of 
Pc_jv2±2i5pc,Arz; withpc^AT^ the percolation threshold ob¬ 
tained from the Newmann-Ziff algorithm, and 5pc,NZ the 
associated error bar. We choose several p values in the 
interval, and by bisecting in between the p’s for which (6) 
becomes linear the earliest, we are able to narrow down 
Pc for precisions as high as six significant digits. For each 
p, we carry out a minimum of 2 x 10^ cluster growths. 
The longer simulations that are carried out for obtaining 
high-precision pc’s can take up to several core-days on a 
computer. 


IV. RESULTS AND DISCUSSION 

We compute the percolation thresholds of the ordered 
coverings of Figs. 1 and 2, as well as the random cov¬ 
erings generated using the algorithms outlined above. 
The findings are summarized in Table I. The Pc’s com¬ 
puted for finite lattices scale linearly with either or 

with L the characteristic length scale of the 



FIG. 7: Application of the epidemic growth method for the 
staggered covering of the triangular lattice. 

lattice, i.e., the number of sites along each dimension. 
This allows us to accurately estimate the percolation 
threshold of an infinite lattice, as depicted in Figs. 1, 2 
and 5. The estimates of Pc on all ordered coverings of the 
square lattice scale linearly with , while all random 

coverings (both on the square and the triangular lattice) 
as well as the two ordered coverings of the triangular 
lattice coverage more rapidly to Pc,oo with the scaling 
For the three coverings with known percola¬ 
tion thresholds, i.e., the parallel and shifted coverings of 
the square lattice (Fig. lb), and the parallel covering of 
the triangular lattice (Fig. 2a), our numerical estimates 
are consistent with the known exact values of Pc, con¬ 
firming the reliability of the Pc’s reported for the other 
ordered and random dimer coverings. 

We use the Pc values obtained from the finite-size scal¬ 
ing analysis as initial guesses for the epidemic growth 
method. As mentioned above, no finite-size correction is 
necessary for the refined Pc’s obtained from the epidemic 
growth method, due to the large lattice sizes considered 
therein. One can therefore test the accuracy of our initial 
finite-size scaling analysis by determining the exponent c 
from: 

\Pc,oo - Pc{L)\ = (7) 

Here Pc,oo and Pc{L) are obtained from the epidemic 
growth algorithm, and the Newman-Ziff algorithm, re¬ 
spectively. Table II summarizes the c values obtained for 
different coverings in this work. These exponents are suf¬ 
ficiently close to the original Yjv and X-VXjv exponents 
used in our initial finite-size scaling analysis. Note that 




















these exponents do not match completely with either of 

or 1 + l/iz because of numerics of considering finite 
number of Pc{Lys in the fit, as well as the statistical un¬ 
certainty of the existing estimates. Also, one cannot rule 
out the possibility of a ‘cross-over’ behavior, which will 
necessitate very large L’s for observing the true asymp¬ 
totic behavior of |pc,oo — PciL)\. 

A. Frustration of the Percolation Transition 

As thoroughly discussed in Section II, dimer covering 
of a lattice leads to a frustration of percolation transi¬ 
tion in the original lattice, and the extent of such frus¬ 
tration can be quantified by computing k = Pc,eq/Pc- 
Since pc is known exactly for the (uncovered) square 
and triangular lattices, k can be readily calculated from 
the computed Pc,eq values. The results are given in Ta¬ 
ble I. The percolation transition is frustrated on dimer 
coverings of both lattices. However, the square lat¬ 
tice tends to be more adversely affected. The paral¬ 
lel and ordered coverings of the square lattice percolate 
at = 0.631966-•• and = 0.626825(2), 

respectively, while the random covering percolates at 
a slightly lower = 0.617713(2). The least af¬ 

fected covering is the shifted covering, which percolates 
at = 0.597296 • • •. Note that the frustration co¬ 

efficient is above 1.19 for all these coverings. Percola¬ 
tion frustration tends to be weaker in the triangular lat¬ 
tice, and K is always below 1.16. Also, the ordered and 
random covering percolation thresholds tend to be much 
closer, and are identical up to the second significant digit. 
This difference can be attributed to higher connectivity 
of the triangular lattice in which the permanently occu¬ 
pied dimer bonds constitute a smaller fraction of total 
bonds. Consequently, the emergence of an infinite clus¬ 
ter is expected to be less sensitive to the a priori occupa¬ 
tion of dimer bonds, as well their specific arrangement. 
It can henceforth be argued that percolation frustration 
will become even weaker in dimer coverings of lattices 
with larger coordination numbers. In two dimensions, 
such lattices can be constructed by adding multiple bonds 
between neighboring sites. Another possibility is to con¬ 
sider three-dimensional lattices, such body-centered cu¬ 
bic and face-centered cubic that have higher coordination 
numbers. 


B. Covering Order and Percolation 

Another interesting question about percolation on 
dimer covered networks is the role of the covering or¬ 
der on percolation. Our findings reveal the difficulty of 
this task, as small changes in the covering can translate 
into significant changes in the percolation threshold. For 
instance, the parallel and shifted coverings of the square 
lattice are almost identical, except for a shift in the dimer 
positions every other row. However, they tend to have 


widely different percolation thresholds. Indeed, their pc’s 
of 0.381966 • • • and 0.347296 • • • tend to fall on the op¬ 
posite ends of the random covering percolation threshold 
^random _ 0.367713(2). This granularity underscores the 
difficulty of identifying an order parameter for character¬ 
izing the role of covering order on percolation, as simple 
order parameters such as (o^) are inadequate for that 
purpose. The origin of this granularity is interesting on 
its own. The shifted covering seems to be a percolation 
‘sweet spot’ on the space of all possible coverings of the 
square lattice. This can be understood by inspecting the 
topologically equivalent networks of ordered and random 
coverings of the square lattice. All these coverings have 
the same coordination number oi z = 6. However, the 
parallel (Fig. la) and the staggered (Fig. lb) covering, 
as well as the random coverings tend to have a large 
number of double bonds in their topologically equivalent 
network, while the shifted covering is exclusively com¬ 
prised of single bonds. Therefore, a single dimer on the 
shifted covering can be connected to a larger number of 
distinct dimers than a similar dimer on the parallel, stag¬ 
gered and random coverings. In other words, the stochas¬ 
tically occupied connections of the shifted network are 
distributed more efficiently as to allow the largest num¬ 
ber of site-to-site connections at a particular occupation 
probability. It can therefore be conjectured that among 
the two-dimensional six-coordinated networks, triangu¬ 
lar networks (whether regular or irregular) will have the 
smallest percolation thresholds. Here, it is necessary 
to mention a related conjecture due to Wierman that 
states that the regular triangular lattice has the highest 
bond percolation threshold among all triangulated net¬ 
works [63]. 

Unlike the square lattice, the difference between or¬ 
dered and random coverings of the triangular lattice is 
very small, and the corresponding percolation thresh¬ 
olds differ only at their third significant digit. Also, the 
two ordered coverings considered here have a higher pc 
than the random covering. In other words, no percola¬ 
tion sweet spots are found among the dimer coverings of 
the triangular lattice. Indeed, no such ordered covering 
might exist as no planar two-dimensional lattice is known 
to have a coordination number of z = 10, which is the 
coordination number of the topological networks of the 
coverings of the triangular lattice. 


C. Comparison with Other Random Networks 

Before comparing the percolation characteristics of 
random dimer coverings with other random networks 
that are defined using the notion of a norm (i.e., distance 
between points in an Euclidean space), it is necessary 
to make a very important distinction. In the norm-based 
random network such as the relative-neighborhood graph 
or the Voronoi tessellation, the coordination number of 
each site is a correlated random variable, although uni¬ 
formity arguments can be made about the average co- 
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TABLE I: Percolation thresholds of dimer-covered lattices. 


Lattice 

Dimer Arrangement 

Pc 

pc,eq 

^ — Pc,eq^lpc 

Square 

Parallel 

0.381966-•• 

0.631966--- 

1.263932- - - 

Square 

Shifted 

0.347296-•• 

0.597296--- 

1.1946- - - 

Square 

Staggered 

0.376825(2) 

0.626825(2) 

1.25365(4) 

Square 

Random 

0.367713(2) 

0.617713(2) 

1.23542(4) 

Triangular 

parallel 

0.237418--- 

0.404085 - - - 

1.163516--- 

Triangular 

Staggered 

0.237497(2) 

0.404162(1) 

1.16374(2) 

Triangular 

Random 

0.235340(1) 

0.402007(1) 

1.15753(2) 


TABLE II: Finite-size scaling exponent, c, for the coverings 
considered in this work. 


Lattice 

Covering 

Exponent 

Square 

Parallel 

0.90 

Square 

Shifted 

0.88 

Square 

Staggered 

0.92 

Square 

Random 

1.97 

Triangular parallel 

1.32 

Triangular Staggered 

1.38 

Triangular Random 

1.72 


ordination numbers based on the stochastic distribution 
of the generating vertices. In the covering-based random 
networks, however, each site can take a finite number of 
coordination numbers that depend on the coordination 
numbers of the original lattice. To be precise, consider 
a dimer bond that connects two sites with coordination 
numbers zi and Z 2 - The corresponding unified dimer 
site on the associated topologically equivalent network 
will have a coordination number of 2:12 = Zi -I- 22 ~ 2- 
This confines the number of allowable coordination num¬ 
bers to a set of known Zij^s. In the square and the 
triangular lattice, all sites have the same 2 , so all sites 
on the covering-based random network will be (2z — 2)- 
coordinated. Therefore, dimer covering can be viewed 
as a systematic way of generating random networks with 
precisely determined coordination numbers. 

Considering the large coordination numbers of the 
covering-based random networks, it is not straightfor¬ 
ward to compare their percolation characteristics with 
those of the known two-dimensional norm-based ran¬ 
dom networks, such as the relative-neighborhood graph, 
Voronoi tessellation and the Gabriel graph that have 2 = 
2.5576 [30], 2 = 3 [26] and 2 = 4 [19], respectively. Not 
surprisingly, these low-coordination graphs tend to have 
higher percolation thresholds, descending in the order of 
increase in 2 . The only comparison can be made between 
the Delaunay triangulation = 0.333069) [26] and the 
random covering of the square lattice Pc = 0.367713(2) 
that are both six-coordinated. Indeed, a random cover¬ 
ing of the square lattice will have a considerable number 
of double bonds, which, as explained above, are likely to 
lead to higher percolation thresholds. In the case of the 


ten-coordinated triangular lattice, no two-dimensional 
random network with a similar coordination number has 
been studied. 


V. CONCLUSION 

In this work, we report the first computational in¬ 
vestigation of bond percolation on ordered and random 
dimer coverings of the square and the triangular lat¬ 
tice. For generating random coverings, we propose a 
novel shuffling algorithm that is conceptually similar to 
an algorithm used for generating proton-disordered ar¬ 
rangements of ice [57]. The percolation transition is 
strongly frustrated in both lattices as a result of dimer 
covering, with the frustration being stronger in the low- 
coordination square lattice. We also observe that there is 
no clear correlation between the existence of long-range 
order in a covering and its percolation characteristics, as 
closely-related ordered coverings can have PcS on oppo¬ 
site ends of the random covering percolation threshold. 
The ordered and random coverings of the triangular lat¬ 
tice considered in this work have very close percolation 
thresholds suggesting that the microscopic detail of the 
covering is less likely to affect the percolation characteris¬ 
tics in highly-coordinated lattices and networks. We also 
compare random coverings with norm-based random lat¬ 
tices such the Voronoi tessellation and Delaunay triangu¬ 
lation both in terms of connectivity and the percolation 
threshold. 

This work paves the way for studying a wide range of 
other interesting questions about percolation on dimer- 
covered lattices. Obviously, it is interesting to study 
bond percolation on dimer coverings of other two- and 
three-dimensional lattices. One interesting question 
there is the sensitivity of percolation frustration to the 
coordination number, dimensionality and the microstruc¬ 
ture of the corresponding lattice. As suggested here, 
and as expected intuitively, frustration is expected to 
be smaller in lattices with higher coordination numbers. 
However, the effect of dimensionality and microstructure 
(i.e., the particular local connectivities in the lattice) is 
less straightforward to predict. Also, the idea of frustra¬ 
tion can be used for understanding and describing the 
behavior of permanently connected systems, such as en- 
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tanglement in polymeric systems. 

Another question of interest is to study dimer coverings 
in which dimer bonds are permanently vacated. Such 
a model can be used for generating random coverings 
with low coordination numbers. For instance, the vacant 
covering of the square and the triangular lattice lattice 
will have a coordination number of z = 3 and z = 5, 
respectively. Finally, one can consider the problem of ri¬ 
mer covering, which is an extension of the dimer covering 
problem. 
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